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In this mainly expository article dealing particularly with recent problems, a computa- 
tional problem of G. Pall related to finite projective geometries is discussed in greater 
detail. Numerical results obtained on SEAC and SWAC are discussed. 



1. Introduction 

All computation on a machine with a fixed number 
of digits is arithmetic modulo some integer and 
therefore can be regarded as number theory. From 
this fact, in "analytic" problems there arise all the 
difficulties associated with "round-off," i.e., neglect 
of insignificant digits. However, there are often 
considerable difficulties in handling actual number 
theoretical problems, e.g., because the numbers may 
grow large rapidly. Codes which permit numbers 
of arbitrary length have been written, but they have 
to be monitored carefully. 

Matrices with integral elements have been studied 
for a very long time and an enormous number of 
problems arise, both theoretical and practical. The 
whole vast classical theory of quadratic forms comes 
under the theory of symmetric matrices; many new 
problems have arisen in recent years. Some of the 
problems can be divided, perhaps, into four topics: 

(1) Inversion (see e.g., J. B. Kosser [1]). 

(2) Eigenvalues and eigenvectors. These play a 
role in the determination of the class number of an 
algebraic number field F by rational methods. For 
it can be shown that the number of ideal classes coin- 
cides with the number of classes of matrices S~*AS, 
where A is a matrix with integral elements which is 
a root of the algebraic equation determining F, and 
S runs through all unimodular matrices with integral 
elements (see O. Taussky [2, 3]). 

(3) Enumeration of matrices with special proper- 
ties. Here we have the various problems concerning 
latin squares, the problems concerning matrices 
with ± 1 as elements and orthogonal rows (see e.g., 
R.E.A.C. Paley [4] and J. Hadamard [5]). The 
study of finite projective geometries has led to new 
problems in this field and some of them will be dis- 
cussed here later. A generalization of these problems 
is given by the study of block designs which are of 
importance in statistical work (see M. Hall [6], 
H. Ryser [7], R. C. Bose [8], W. S. Connor and 
W. H. Clatworthy [9]). 

(4) Maximization problems. These come actually 
also under the problems discussed in (3), but have 
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been studied rather extensively, particularly in recent 
years. A prototype among these is the assignment 
problem: given an nXn matrix (a tj ), determine the 
permutation P(i) for which 

is maximum (see e.g., T. S. Motzkin [11], C. B. 
Tompkins [12], H. W. Kuhn and A. W. Tucker [13]). 

2. Some Problems 

However, many other problems arise which cannot 
easily be brought under these headings. Of these 
we mention a few. Matrices with rational integral 
coefficients are a generalization of rational integers 
and any number theoretical process can be applied 
to them: e.g., any two nXn matrices of this type 
have a greatest common left divisor, a greatest 
common right divisor, a least common left multiple, 
a least common right multiple. If D is the greatest 
common left divisor of A,B, then 

D=AP+BQ 

can be solved where P, Q are again nXn matrices 
with rational integral coefficients. Diophantine 
equations can be studied, even the Fermat equation 

x n -\-y n =z n 

can be considered for matrices. 

Completely new problems arise as well because of 
the noncommutativity of matrices. A quite elemen- 

(a b\ 
tary example is to express a matrix I , with a, 

\c dj 
b, c, d rational integers and ad—bc= 1 in terms of the 

/l 1\ /0 -1\ 

generators and . A much more dif- 

\0 1/ \1 0/ 

ncult problem is to express (when possible) a positive 
definite nXn matrix C of rational integers in the 
form C=XX' where X is again an nXn matrix of 
rational integers and X f is its transpose. If the 
determinant |C| = 1 and n<8 such a decomposition 
is known to be possible always. Actually to find 
X can be a difficult computational problem (see [7]). 
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3. Hermite Normal Forms 

For many of the problems concerning matrices of 
rational integers the so-called Hermite normal form of 
a matrix A plays an important role. It is a triangular 
matrix obtained by multiplying A with a suitable 
unimodular matrix U of integers and is of the form. 



hn h\ 2 . . hi r 

h 22 . . h? Tl 
UA= 



L.o o 



o K 



where h u >0 and 0<h ik <^h u (i<^k). 

Another important normal form associated with 
A is the Smith normal form which is a diagonal 
matrix. It is obtained by multiplying A with two 
unimodular matrices U, V so that 



UAV= 



Oj 



and h 1 \h 2 



\h P 



4. Application of the Hermite Normal Form 
to Finite Projective Planes 

An important use of the Hermite normal form was 
made in the study of finite projective plane geometries 
(FPP). An FPP is a set of k elements called points 
which are distributed over k lines such that exactly 
ft+1 points lie on each line and exactly n+1 lines 
pass through a given point. Two distinct points 
are contained in one and only one line; two distinct 
lines contain one and only one point in common. We 
have k=n 2J \-n+l. The number n is called the 
order of the geometry. With such a geometry there 
is associated the so-called kXk incidence matrix 



A= {a i3 ) 



where 



line. 



a tj = 1 if the ith point lies on the jth line, 
aij=0 if the ith point does not lie on the jth 



Hence any two rows (columns) have a one in common 
in exactly one column (row) and there are n-\-l 
elements a u =l in each row (column). 

The incidence matrix A is a matrix of integers. 
If the lines (points) are numbered in a different order 
then A is replaced by PA (AQ) where P and Q are 
permutation matrices. The following two properties 
of A are particularly important: 

det A=±n< n2 + nU2 (n-{-l), 
AA'=A'A=nI+J, 

where / is the unit matrix and J the matrix consisting 
of l's only. 

G. Pall studied the Smith and Hermite normal 
forms of A (see [14]). He showed that for square 
free values of n the diagonal elements in the Smith 
form are 

1, 1, . . . , 1, n, . . . , n } n(n-\-l) 

and the diagonal elements of the Hermite normal 
form are 1, or divisors of n, or (n-^-l) times a divisor 
of n. 

Two different geometries cannot have the same 
Hermite matrix. Thus the Hermite matrix of a 
geometry is rather important. G. Pall further 
investigated the ^-rank of A, i.e., the maximum 
number of linearly independent rows of A, mod p, 
where p is a prime factor of n. If p 2 is not a factor 
of n (e.g., if n is square free) then the ^-rank of A 
is exactly 



1 



(n 2 +n) + l- 



For all other p 



p-r&nkK- (n 2 -\-ri)-\-\ 

L 

holds. As special examples the 2-rank for the 
(unique) geometry corresponding to n=4: was com- 
puted. It is 10. The 3-rank corresponding to 
n=9 presented a more difficult task. There are 
two different geometries of order 9. The 91X91 
incidence matrices of these geometries have as 
3-rank 37 and 41 while ^(r^ 2 +^) + l=46. These 
values of the 3-rank were computed on SEAC and 
on SWAC. The computation on SEAC was carried 
out by M. Newman; it took 20 hr. The Hermite 
normal forms of the geometries of order n = 5 and 
n=7 were also computed; it took 3 hr for this task. 

5. Further Facts About the Incidence 
Matrix 3 

Incidence matrices have been used to study the 
still unsolved problem: For what value of n is there 
a finite projective plane? So far it is known that 
for n=p r , p a prime number, there is always a 
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finite projective plane. This plane is uniquely 
determined for n = 2, 3, 4, 5, 7, 8. This was estab- 
lished for n=7 by M. Hall [16] and for n=S by 
computations carried out on SWAC by M. Hall, 
J. D. Swift and R. J. Walker [17]. 

For n = 6 there is no FPP whatsoever (see R. H. 
Brack and H. J. Ryser [18], S. Chowla and H. J. 
ftyser [19]). On the other hand for n=9 it is known 
that there exist two different geometries; one of these 
satisfies the axiom of Desargues, the other does not. 
For n=10 it is not known whether a FPP exists. 

A valuable contribution towards the solution of 
the problem of determining the possible orders of a 
FPP is given by the following theorem of Ryser [20]: 

Let n be odd and let k=n 2J r n J r l. Assume nl+J 
= XX' where X is a matrix oj integers. Then there 
exists a FPP with n + 1 points on each line. 

However, as pointed out earlier, it is a rather 
difficult computational problem of deciding whether 
a given positive definite matrix of high order splits 
up in the form XX* with X an integral matrix. 
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